Photonic excess noise and wave localization 
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This is a theory for the effect of locahzation on the super-Poissonian noise of radiation propagating 
through an absorbing disordered waveguide. Locahzation suppresses both the mean photon current I 
and the noise power P, but the Fano factor P/I is found to remain unaffected. For strong absorption 
the Fano factor has the universal value 1 + §/ (with / the Bose- Einstein function), regardless of 
whether the waveguide is long or short compared to the localization length. 

PACS numbers; 42.25.Dd, 42.50.Ar 
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The coherent radiation from a laser has Poisson statis- 
tics . Its noise power Ppoisson equals the mean current 
/ (in units of photons per second) . Elastic scattering has 
no effect on the noise, because the radiation remains in 
a coherent state. The coherent state is degraded by ab- 
sorption, resulting in an excess noise P — Ppoisson > 
[^. The Fano factor P/Ppoisson deviates from unity by 
an amount proportional to the Bose-Einstein function /. 
It is a small effect (/ ~ 10~^ at room temperature for 
infrared frequencies), but of interest because of its fun- 
damental origin: The excess noise is required to preserve 
the canonical commutation relations of the electromag- 
netic field in an absorbing dielectric |^-^ . 

The interference of multiply scattered waves may lead 
to localization [Q . Localization suppresses both the mean 
current and the fluctuations — on top of the suppression 
due to absorption. Localization is readily observed in a 
waveguide geometry , where it sets in once the length 
L of the waveguide becomes longer than the localization 
length ^ ~ Nl (with I the mean free path and N the num- 
ber of propagating modes). Typically, £, is much larger 
than the absorption length ^a, so that localization and 
absorption coexist. The interplay of absorption and lo- 
calization has been studied previously for the mean cur- 
rent Here we go beyond these studies to include 
the current fluctuations. 

It is instructive to contrast the super-Poissonian pho- 
tonic noise with the sub-Poissonian electronic analogue. 
In the case of electrical conduction through a disordered 
wire, the (zero-temperature) noise power is smaller than 
the Poisson value as a result of Fermi-Dirac statistics. 
The reduction is a factor 1 /3 in the absence of localiza- 
tion (l^Jl^ . The effect of localization is to restore Pois- 
son statistics, so that the Fano factor increases from 1/3 
to 1 when L becomes larger than ^. What we will show 
in this paper is that the photonic excess noise responds 
entirely differently to localization: Although localization 
suppresses P and /, the Fano factor remains unaffected, 
equal to t he valu e 1 -I- |/ obtained in the absence of lo- 
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tion of the problem. The noise power 
P= / dt6I{0)SI{t) 



(1) 



quantifies the size of the time-dependent fluctuations 
of the photon current I{t) — I + 61{t). (The bar 
indicates an average over many measurements on 
the same system.) For a Poisson process, the power 
Ppoisson ~ I equals the mean current and the Fano fac- 
tor J- = P/Ppoisson equals unity. We consider monochro- 
matic radiation (frequency luq) incident in a single mode 
(labeled toq) on a waveguide containing a disordered 
medium (at temperature T). (See Fig. |l|.) The incident 
radiation has Fano factor J^n. We wish to know how the 
Fano factor changes as the radiation propagates through 
the waveguide. 




FIG. 1. Monochromatic radiation (thick arrow) is inci- 
dent on a disordered absorbing medium (shaded), embedded 
in a waveguide. The transmitted radiation is measured by a 
photodetector. 

Starting point of our investigation is a formula that 
relates the Fano factor to the scattering matrix of the 
medium |15|, 
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Let us begin our analysis with a more precise formula- 



(We have assumed detection with quantum efficiency 1 
in a narrow frequency interval around loq.) The function 
f{u!,T) = [exp{hu!/kT) — is the Bose-Einstein func- 
tion. The transmission matrix t and the reflection matrix 
r are N x N matrices, with N the number of propagating 
modes at frequency ujo. The term proportional to / in 
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Eq. (g) is the excess noise. For a unitary scattering ma- 
trix, rr^ +tt^ equals the unit matrix 1, hence the excess 
noise vanishes. 

In what foUows we will assume that the incident radi- 
ation is in a coherent state, so that J^m = 1 and the devi- 
ation of J-' from unity is entirely due to the excess noise. 
Since the Bose-Einstein function at room temperature is 
negligibly small at optical frequencies, one would need to 
use the coherent radiation from an infrared or microwave 
laser. Alternatively, one could use a non-coherent source 
and extract the excess noise contribution by subtracting 
the noise at low temperature from that at room temper- 
ature. 

The absorbing disordered waveguide is characterized 
by four length scales: the wavelength A, the mean free 
path for scattering I, the absorption length ^a, and the 
localization length ^ = {N+l)l. We assume the ordering 
of length scales A <C / <C <C which is the usual situ- 
ation 1^. We ask for the average (J-) of the Fano factor, 
averaged over an ensemble of waveguides with different 
realizations of the disorder. For L ^ ^a we may neglect 
the matrix tt'' with respect to 1 in Eq. (||), so that the 
expression for the Fano factor (with JFj„ = 1) takes the 
form 

.F=l + 2/(l-Ci), C,^ [^^^"ry^^"'""" . (3) 

In the absence of localization, for L ^ ^, one can sim- 
plify the calculation of (JF) by averaging separately the 
numerator and denominator in the coefficient Ci, since 
the sample-to-sample fluctuations are small. This diffu- 
sive regime was studied in Refs. p^ , p^ . Such a simpli- 
fication is no longer possible in the localized regime and 
we should proceed differently. 

We follow the general approach of Ref. ||l^, by con- 
sidering the change in upon attaching a short segment 
of length 6L to one end of the waveguide. Transmission 
and reflection matrices are changed to leading order in 
SL according to 

t tsL{l + rrsL)t, r^r'gL+tsL{l + rrsL)rtJi^, (4) 

where the superscript T indicates the transpose of a ma- 
trix. (Because of reciprocity the transmission matrix 
from left to right equals the transpose of the transmission 
matrix from right to left.) The transmission and reflec- 
tion matrices tsL, fsLi of the short segment have zero 
mean and variances 

kl viL\ kl FSLi mn I 

= (4m<^/n + 5kn5lm)5L/i, (5a) 
{[tMULY^n) = N~Hk,n5ln{l " 5L/1 - SL/Q, (5b) 

where la = is the ballistic absorption length. All 

covariances vanish. Substituting Eq. (|^) into Eq. (||) and 
averaging we find the evolution equation 
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where we have defined pp = tr (1 — rr'^ Y . 

For L ^a we may replace the average of the prod- 
uct (Cipi) by the product of averages (Ci)(pi), because 
|l2| statistical correlations with traces that involve reflec- 
tion matrices only are of relative order ^a/C — which we 
have assumed to be <C 1. The moments of the reflection 
matrix are given for i ^ ^a by 
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hence they are ^ 1 and also ^ SI . We may therefore 
neglect the terms in the second, third, and fourth line of 
Eq. (0). What remains is the differential equation 



d{Ci) 
dL 



-2(Ci)(pi) + (p2), 
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which for L ^ has the solution 

(P2) 
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We conclude that the average Fano factor {!F) = 1 + 
2/(1 - (Cl)) H- §/ for L > ^a, regardless of whether 
L is small or large compared to ^. 




FIG. 2. Length dependence of the average Fano factor, 
computed from Eq. ^ with Tin ~ 1. The data points re- 
sult from a numerical simulation for an absorbing disordered 
waveguide with A'^ = 50 propagating modes. Arrows indicate 
the absorption length and the localization length ^. The 
average Fano factor is not affected by localization. 
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FIG. 3. Dependence of the average and variance of the 
Fano factor on the number A'^ of propagating modes, for fixed 
length L — 260 Z = 38.5 of the waveguide. The length is 
larger than the localization length ^ — {N + 1)1 for all data 
points. The dashed lines extrapolate to the theoretical expec- 
tation for 1/iV ^ 0. 

To support this analytical calculation we have car- 
ried out numerical simulations. The absorbing disor- 
dered waveguide is modeled by a two-dimensional square 
lattice (lattice constant a). The dielectric constant e 
has a real part that fluctuates from site to site and a 
non- fluctuating imaginary part. The multiple scatter- 
ing of a scalar wave 5* is described by discretizing the 
Helmholtz equation [V^ + ((jJo/c)^£]^' = and comput- 
ing the transmission and reflection matrices using the re- 
cursive Green function technique |p^ . The mean free 
path I — 20 a and the absorption length = 135 a 
are determined from the average transmission probability 
N-'^{tTtt^) = Z/^aSinh(L/^a) in the diffusive regime [jl2l. 
Averages were performed over the N/2 modes mo near 
normal incidence and over some 10^ — 10^ realizations of 
the disorder. Results are shown in Figs. || and ||. 

The length dependence of the average Fano factor is 
plotted in Fig. 2, for = 50 and L ranging from to 
2^. Clearly, localization has no effect. The limiting value 
of f~^{^ — 1) resulting from this simulation is slightly 
smaller than the value 3/2 predicted by the analytical 
theory for N 1. The iV-dependence of {T) in the lo- 
calized regime is shown in Fig. 3. A line through the 
data points extrapolates to the theoretical expectation 
- 1> ^ 3/2 for N ^ oo. Fig. 3 also shows the 
variance of the Fano factor. The variance extrapolates 
to for oo, indicating that T = P/I becomes 

self-averaging for large N . This is in contrast to P and 
/ themselves, which fluctuate strongly in the localized 



regime. 

In conclusion, we have demonstrated that localization 
of radiation in an absorbing disordered waveguide has 
no effect on the ratio of the excess noise and the mean 
current. In the limit of a large number of propagating 
modes, this ratio is self-averaging and takes on the uni- 
versal value of 3/2 times the Bose-Einstein function. Ob- 
servation of this photonic analogue of the universal 1/3 
reduction of electronic shot noise presents an experimen- 
tal challenge. 
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